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The matrix element for the decay of orthopositronium to three photons can be expressed in 
terms of three independent amplitudes. We describe the analytic evaluation of these ampli- 
tudes, both to lowest order and with the inclusion of all one-loop corrections. We use these 
amplitudes to find precise values for the one-loop correction to the orthopositronium decay rate 
Ti = — 10. 286606(10) (a/7r)rLO, and for the order-a 2 "square" correction to the decay rate 
^(square) = 28.860(2)(a/7r) 2 rLO, where Flo is the lowest order rate. We give in explicit form 
the function describing the one-loop correction to the distribution in phase space of the final state 
■ photons. 
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£ ' I. INTRODUCTION 

3 ' 

Positronium, the electron-positron bound state, is well-suited for probing many fundamental aspects of particle 
• physics. The physics of positronium is governed almost exclusively by the electromagnetic force — weak interaction 
effects are negligible compared to present experimental and theoretical uncertainties 0, 0, 0, ■ As a consequence, 
positronium is an ideal system for testing QED through high precision comparison between experimental and theo- 
retical results for energy levels and decay rates. The states of positronium are eigenstates of the charge conjugation 
and parity operators C and P, so positronium can be used to test the discrete symmetries C, P, and T and combina- 
tions thereof. Positronium has been the focus of many past and ongoing attempts to observe physics beyond the 
standard model. 0,0 In this work we focus on the decay of spin-1 orthopositronium to three photons. 

The orthopositronium decay rate has been the subject of continuing experimental and theoretical work since the 
ly-^ first measurement by Deutsch in 1951. 0] A summary of all experimental and theoretical results has been given by 
Adkins, Fell, and Sapirstein ^} an( l updated with commentary by Rubbia Q and by Sillou [ll|. By 1990 it was 
apparent that there was an "orthopositronium lifetime puzzle" , as the most precise experimental determinations (gas 
UJ and vacuum 0] results from the Michigan group) were in disagreement with theory 0, by several standard 
deviations. Many experiments were mounted to look for exotic decays of orthopositronium in an attempt to resolve 
Qh, the discrepancy, @, 0, 0, 0, Newer, somewhat less precise powder results from the Tokyo group in 1995 [2(j 
and 2000 |2l| were consistent with theory and inconsistent with the earlier Michigan results. In 2000 the calculation 
of all Q(ce 2 ) corrections to the decay rate were c omp leted. Including yet higher order logarithmic corrections 

as well, [22, 0, E3 the theoretical prediction is |l0j 

r(theory) = 7.039979(ll)^s" 1 . (1) 

The 0(a 2 ) correction was found to be not unusually large, leaving the discrepancy with the Michigan results intact. 
Finally, in 2003 the lifetime puzzle was resolved by two new high-precision results from the Tokyo [2(| and Michigan 
[27| groups: 

T(Tokyo) = 7.0396(12 stat.)(ll syst.)^ 1 (2a) 
F(Michigan) = 7.0404(10 stat.) (8 syst.)^ 1 , (2b) 

consistent with each other and with theory. 

The resolution of the o-Ps lifetime puzzle does not decrease the long-term usefulness of positronium decay as 
a probe of fundamental physics. Ongoing and proposed experiments involving positronium decay include those of 
Refs. [a, H H3i HE • One challenge is to improve the experimental precision of the o-Ps decay rate (currently about 
200 ppm) to a level closer to the present theoretical value (about 2 ppm). The 0{a 2 ) contribution to that rate is 250 
ppm, so improved experimental precision will be required in order to test the 0(a 2 ) calculated result. 

In this work we describe an analytic evaluation of the one- loop o-Ps— > 37 decay amplitudes. We use these amplitudes 
to obtain a precise value for the 0(a) decay rate contribution, and also to calculate the part of the 0(a 2 ) correction 
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coming from the square of the one-loop amplitudes. These results have been reported already |3fJ — here we give 
further details. We also supply an explicit analytic expression for the 0(a) differential decay rate in terms of photon 
energy variables. From the differential decay rate it is easy to obtain the 0(a) corrected one-photon energy spectrum. 
(This energy spectrum, calculated more laboriously by numerical methods, has been useful in developing simulations 
of experimental arrangements. |20ll26| ) 

We ad apt the formalism of covariant decay amplitudes, originally developed for the study of Z boson decay to three 
photons, |31j to the case of o-Ps— > 37. In Sec. [H] we use the extensive symmetries of the decay tensor to show that 
there are only three independent amplitudes for the o-Ps— > 37 decay. In Sec. IIIII we express the decay amplitudes in 
terms of helicity variables since the spin sums are most convenient in this form. In Sec. II VI the integral for the decay 
rate is reduced to its minimal two-dimensional form. In Sec.lVlt he p receeding formalism is applied to the lowest-order 
decay process and the lowest-order decay rate of Ore and Powell [33] is reproduced. In Sec. I VII the method of Passarino 
and Vcltman 33] for evaluating one-loop integrals is developed. In Sec. IVIII the one-loop calculation is described. 
Finally, in Sec. I Villi our results for the 0(a) and part of the 0(a 2 ) decay rates are given. The Appendix contains our 
explicit form for the one-loop decay distribution. 



II. SYMMETRIES OF THE DECAY TENSOR 

The decay of the massive vector particle orthopositronium to three photons is described by the matrix element 



M = e* 1 ^e* 2fl2 e* 3 ^e a M^°<*(k 1 ,k 2 ,k 3 ) , (3) 

where the three photons have momenta fcj and polarizations ej, and the positronium atom has momentum P = 
k\ + k 2 + k 3 and polarization e. The decay tensor is a linear combination of terms like k^ 1 k£ 2 fc£' 3 k%, k^ 1 k^g^ 3 " 1 , and 
gMiM2gM3a xhe most general such tensor has 81 terms of the first type, 54 of the second, and 3 of the third. However, 
gauge invariance and Bose symmetry reduce the number of independent contributions to only three |3lT ] . We review 
the argument below. 

Because the decay tensor is always contracted with physical polarization vectors of on-shcll photons, which satisfy 
Ea^a = (f° r a = 1 , 2, or 3), we can drop terms containing factors of fc^ 1 , k 2 2 , and k^ 3 . This leaves only 24 terms 
of the first type, 30 of the second, and still 3 of the third. 

By Bose symmetry, the tensor M is totally symmetric under photon interchange. This means, for the interchange 
of photons 1 and 2, that 

M^^ a (k 1 ,k 2 ,k 3 )=M tl2 ^ sa (k2,h,k 3 ) ■ (4) 

This symmetry leaves only four independent terms of the first type, six of the second, and one of the third. The decay 
tensor can be written in the manifestly symmetric way 

M»^ a (k 1 ,k 2 ,k 3 )=Y,M msa (ki,k 2 ,k 3 ) , (5) 
s 3 

where the sum is over the six photon permutations, and the tensor M. has the form 

M f * 1 ^ sa (ki,k 2 ,k 3 ) = ai(fci,A: 2 ,fc 3 )fcf fcf fcf kf + a 2 {k u k 2 ,k 3 )k% 1 k^k^k« 
+ a 3 (fci, k 2 , k 3 )k^ fcffcf fcf + a 4 (fci, k 2 , fc 3 )£;ffcf fcffc? 

+ 63(fci,fc 2 ,fc 3 )fc 3 11 W 2 ^ +6 4 (fci,fc2,fc 3 )fcf fc£V lQ 
+ h(k 1: k 2 , k 3 )kfk^g^ a + b 6 (k u k 2 , k 3 )k^k^g^ a 

+ c(h,k 2 ,k 3 )g^ a g^ . (6) 

The quantities a,i(ki, k 2 , k 3 ), bi{k\, k 2 , k 3 ), and c(fci, k 2 , k 3 ) are scalar functions of their arguments, and 65, be, and c 
are symmetric under the interchange k 2 <-> k 3 . 

Gauge invariance requires that the tensor M be transverse 

k^ x M^^ a (k u k 2) k 3 ) = , (7) 

with similar relations holding for contractions with k 2fl2 and k 3fl3 . The condition of Eq. (JJJ) provides 13 independent 
relations among the 19 variables a^fci, k 2 , k 3 ), aj(fci, fe 3 , k 2 ), bi(ki,k 2 ,k 3 ), b t (ki, k 3 , k 2 ), b 5 (ki,k 2 ,k 3 ), b 6 (ki,k 2 ,k 3 ), 
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and c(ki,k 2 ,k 3 ) (with i = 1,2,3,4), which lead (using permutation symmetry) to three independent solutions. So, 
the tensor M. can be expressed in terms of three scalar functions A\, A 2 , and A 3 as 

1 /k^k^ 3 \ / k^ 2 k^ 2 \ 

1 / k"^ 1 \ / k^ 2 k^ 3 \ 

+ *- fe '«ir5(^-'")(li-»~) ■ <» 

The amplitudes A\, A 2 , and A 3 can be identified by writing the decay tensor as 

M"™ a (fci, k 2 , fcs) = - i4i(fci, fe, k 3 )k^k^ 2 fcffc? [(fci • fc 3 ) 2 (fc! • fea)] _1 

+ A 2 {k 1 ,k2,k 3 )k^ 1 k^ fcf fcf [(fci-feKfe-fcsK^-*!)]" 1 

+ A 3 (fc 1 ,A :2 ,fc3)fc3 !lfc 3 2fc 2 3fc ?[( fc l ' fc 3) 2 (fc 2 • fc 3 )] _1 

+ ... • (9) 
One finds Ai, A 2 , and A 3 by taking the coefficients of fcffcf fcf jfef , fcffcf fcf fcf , and fc^/cffcf fcf . 

III. HELICITY AMPLITUDES 

The formula for the o-Ps — > 37 decay rate involves the absolute square of the decay matrix element summed over 
final state spins and averaged over the initial state spin: 

This is a Lorentz invariant quantity, and can be calculated in any frame. It is convenient to calculate it in a two-photon 
rest frame. 

Since we will use the orthopositronium center of mass frame for the decay rate integration, it is useful to express 
our results in terms of invariant variables. A convenient set is given by the Mandelstam variables, which are defined 

by 

Sij = Sji = [ki -\- kj*) = 2ki • kj , (1 1) 

and satisfy 

Sij + s jk + s kl = Mp s , (12) 
where Mp s here is the orthopositronium mass and {i, j, k} is any permutation of {1, 2, 3}. Bar variables are defined 

by 

= Mp s - = s ik + s jk . (13) 

They satisfy 

% + Sj k + s kl = 2M| S . (14) 

We note that each Sjj and s^ is non-negative. 

We calculate \M\ 2 in the k\k 2 rest frame. The photon and positronium momentum vectors in (E,p x ,p y ,p z ) notation 
are given by 

fci = (fc,0,0,fc) 

k 2 = (fc,0,0,-fc) 

&3 = {q,qsin9,0,qcos9) 



P = (E,qsm6,0,qcos6) , 
where E 2 = q 2 + Mp s . The k\k 2 rest frame kinematic variables are given in terms of 



k 
Q 

E 
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The helicity vectors for photon 1 are 



et = ^(0,-1, -i,0) 
er = ^(0,1, -i,0) 

For photon 2 we rotate these by 180° around the y axis using 

-10 
R 2 = | 1 
-1. 



to find 



For photon 3 we rotate using 



R 2 ef 



cos 9 sin 
H, = | 10 

•sin # cos , 



to find 



eg = R^e^ = — ^=(0, — cos#, — i, sm9) , 
1 

eg = i?3e^ = — p(0, cos#, — i, — sm8) . 
v2 

The positronium spin ±1 helicity vectors are the same as those for photon 3: 

c Ps " e 3 ' 

while the positronium spin helicity vector is 

= -^-(q,E sine, 0,Ecos6) . 

The helicity amplitudes are defined by 

^AiA 2 A 3 A — e lf 11 e 2^2 e 3 f i 3 e Psa m [Kl,K 2 ,K 3 ) 

There are nine independent helicity amplitudes with Ai = +. They are 

r-Ai(123) + Ai(132) , A 2 {123) - A 2 (132) 
M ++++ - l \ = 1 = 

L «12 Sl2 

S23 Ag(132) Ag(312) | 2) -| 

S12S12 S 23 / 

= ^ ^(123)^^(132) | -A 2 (123) + ^ 2 (132) 

I S\ 2 S12 



■5 



Sl3 S12S12 



(l<-2)} , (25b) 
M ++ _+ = 2 {-^ffi>+(l~2)} , (25c) 

Af ++ __ = 21^2) + -A 2 (123)-A 2 (132) _ *(312) + _ | 

L S12 S23 s 23 J 

rAi(123) -^2(123) - A 2 (132) -A 2 (312)- A 2 (321) 
2s - -|- - -(- 

L S12 S12 S12 

_ ^(213) £23^3(231) I 

S23 S12S12 J 

r ^(123) | 8 i 3 (-^2(123)-A 2 (132)) si 3 A 3 (231) | 

t S12 S 23 Si2 S12S12 J 

2A ^3(^(123) -A!(132)) _ r 13 (-A 2 (123) + A 2 (132)) 
MS"! I7o ^1(312) + 



M + _ + 



+++0 



M£ s I S12 S12 

+5 12 A 2 (312) + S12S2 3A 3 (123) + (Mg s + .12)^3^23^3(132) 

Sl3 Si 2 Si 2 

+ gl2gl 3A 3 (312) _ (i ^ 2) | (25g) 



S23 



M ++ _ = ^( ^(132) - r 23 A l( 231) + + Mm _ ^ (231)) 



_£i2£] 

S23 Sl3 

M + . +0 = ^-{--ffl_ 5l2Al(321) 
Mp s S12 

+ 



^12^13^3(312) ^ si 2 s 23 A 3 (321) J (25h) 



(Mg s + si 2 )si3(^ 2 (123) + A 2 (132)) + MMm) + Aa(32i)) 



S12 

S23 S i2 Si 2 



^12^13^3(213) | (Mg s + si 2 )si 3 s 23 A 3 (231) j ^ 



where we have used the abbreviated notation Ai(abc) = A i (k a7 ki )7 k c ) and the definitions 

Tij — Afp s Sjj Sik s jk 1 



A = L Mk . (26) 
V 2si 2 si 3 s 23 

The other three Ai = + amplitudes are related to the previous ones by 

M+—+ = M+_+_(1^2) , 

M+___ = M+_++(l <- 2) , 

M+__ = M+_ +0 (l <-» 2) . (27) 
The Ai = — amplitudes are given by the parity relations 

M-\ 2 \ 3 ± = M + _a 2 -A 3 t > 

M_ A2A3 o = -M+_ A2 _ A3 o . (28) 
The squared decay matrix element can be written as 

W?=\ E l M +w| 2 • (29) 

IV. THE DECAY RATE INTEGRAL 

We will calculate the o-Ps — > 37 decay rate integral in the positronium center of mass frame. The decay rate 
integral is given by 



1 1 f d z k x d 3 k 2 d 3 k 



r = 312^7 (^(^(^^ P - fcl - fc2 -^' M ' 2 ^ 



6 



\AAAAf 



vAAAAj 



vAA/W- 



FIG. 1: The lowest order o-Ps — > 37 decay graph. The factor on the right represents the initial spin-1 wave function. 



where Wi = \ki\ are the photon energies. Of the nine variables in k\, k 2 , k 3 , four are determined in terms of the others 
by energy-momentum conservation 



U!l + LU 2 + UJ3 = MPs 

ki + k 2 + k 3 = . 



(31) 



Three variables describe the orientation in space of the decay plane. The remaining two variables describe the relative 
orientation of the photons in the decay plane. We will use the energies of two of the photons for this last pair 
of variables. Each photon can have any energy between and W = Mp s /2. We find it convenient to introduce 
dimensionless variables Xi = LUi/W which satisfy < x\ < 1, X\ -\-x% + x$ = 2 and are given in terms of invariants by 
Xi — Sjk/Mp s . In terms of the x's, one has 



W 



768Wo d n^ 2|M|2 



(32) 



V. THE LOWEST ORDER DECAY RATE 



The lowest order decay amplitude is given by 



s 3 



-(-• ie 7 e|)- 



-(-iejel)^] , 



3 ' 1 (-P/2 + k 3 )-m y ^'^' 1 (P/2-k 1 )-m' 
where the sum is over the six permutations of the final state photons. The wave function factor is given by 

which contains the spin-1 spin factor, a normalization factor, and the wave function at contact 



8tt 



We write 



V) =\(lN + l)7e(7N-l) 
for the positronium spin factor where N = P/(2W). The lowest order decay amplitude (see Fig. 1) becomes 



J2 ^tr[ 7 e*(- 7 i? 3 + 1)7^(7^1 + 1M(7^ + l) 7 e( 7 iV - 1)] 



-lira 



4 e — ' Xi x 3 

s 3 1 d 



where Ri = N — Ki, Ki = ki/W and W ' « m. The lowest order decay tensor has the corresponding form 



JW LO 



V — !— trlj^i-jRs + l)y»(7iZi + 1)7^ ( 7 iV + 1) 7 Q (7^V - 1)1 



4 " £1X3 
03 



(33) 



(34) 



(35) 



(36) 



(37) 



(38) 
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We replace N by (Ki + K2 + K^)/2 and expand this out, and identify the Ai functions by use of Eq. JSJ). The lowest 
order functions are 

A\ (x 1 ,x 2 ,x 3 ) = , 

AhOf \ if 2 3 X \X2X 3 

X1X 2 X 3 

A^°( Xl ,x 2 ,x 3 ) = , (39) 
where Xi — Sjk/Mp s and x% — 1 — x% = Sjk/Mp s . Clearly, A\° is a factor in each hclicity amplitude. One has 

1U +-++ 



lu +-+o 



-X3 






X1X2 




m 


• 


M™ 




-1 1 


x 3 x 3 ; 


-x 2 


A\° 




X1X3 


m 2 


• 


-Xl 


A%° 




X2X3 


m 2 


1 


X2 I 




A\° 




X1X2X3 


m 2 


Xl 1 




A\° 


~ 3 y 


X1X2X3 


m 2 



x 3 V x lX2 x 3 

with all other Mffi A A amplitudes equal to zero. One finds that 

^^A 6 {(iL) 2 + (A) 2 + (A) 2 } . (41) 

6 l\X2X 3 / \XlX 3 / \XlX2' > 

The lowest order decay rate is the Ore and Powell result [13 

r LO = ^ 6 f 1 dx 1 f ^{(A) 2 + (A) 2 + (iL) 2 } 

9ir J Q Ji-xi t\ X2X3 / \xix 3 J \ Xl x 2 / > 
= l(^_9)mo 6 . (42) 

97T 



VI. ONE-LOOP INTEGRALS 

We used the method of Passarino and Veltman to evaluate the one-loop integrals. Since this approach is widely 
used, and lengthy to describe in detail, we will just list the one-loop integrals that are required but only work the scalar 
integrals out in detail. The Passarino- Veltman method will be illustrated in the case of the three-point functions. 

The general definition of the one-loop form factors is through 

{X^X^X^,--} = n 2e J \dq)n{l,qv.,q»q v ,- • } 

x [(~q 2 + m 2 1 )(-(q+p 1 ) 2 + m 2 2)(-(q+p 1 +p 2 ) 2 + ml)---y 1 . (43) 

Ultraviolet divergences are controlled through dimensional regularization with n = 4 — 2e the dimensionality of 
spacetime. We define (dq)',[ = d n qj ii-K n l 2 ). The quantity fj, is a reference mass introduced with the regularization 
which we take to be equal to the electron mass m. Functional dependences on the masses and momenta are indicated 
by X(m 1 ,m 2 ,m 3 , ■■■ ;p 1 ,p 2 , ■ ■ ■)■ 

The one-point function is trivially evaluated: 

A(mi) = m 2£ / (rig)" 2 \ 

J (~q 2 + mf) 

= -m\Y{e) + A{mx) + 0(e) (44) 

where 



A(m\) = —to 2 [l — ln(TO 2 /m 2 )] 



(45) 
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The two-point functions are defined by 

{B ,B^,B^} = m 2e I \dq)'^{l,q^q^}[(-q 2 + m 2 ) (- (q + p) 2 + m^f ] . (46) 



The scalar function Bq is 

= r{e) + B (m 1 ,m 2 ;p) (47) 

where A = (f — x)m\ + xm 2 , — x(l — x)p 2 and 

B (m 1 ,m 2 ;p) = -J dx\n(A/m 2 ) . (48) 
All parametric integrals will be taken between the limits and f . The cases of interest are 

B (0,m;p) = 2+^— £la(l-p) (49) 

where p = p 2 /m 2 , and 

B (m,m;p) =2|l- ~ arctan } > ( 50 ) 

valid for < p < 4. 

The three-point functions are defined by 

{<7 , (7^,(7^,(7^} = m 2e J(dq)" { 1 , g M , , ^q^Qq } 

x [(- g 2 +m?)(-( g +p 1 ) 2 + m|)(-(g + p 1 +p 2 ) 2 + m2)] _1 . (51) 
The general forms for C^, C^ u , and are 

<7 M = Pi M <7n +p 2AI Ci2 , (52a) 

Cpv = V\ i iV\vC2\ +P2fiP2uC 2 2 + {PlP 2 }^C 23 + g^ v C 2 A , (52b) 
C^a = Pl^Plt/PlaCai + P2 t iP2uP2aC 3 2 + {PlPl^l/ii/aC^ 

+{j3iP2P2}^ 1 /q(7 3 4 + {pi5j> ra C35 + {p 2 g}^ a C 36 , (52c) 

where 

{pk}^ = p^k u + k^p 
{ppk} 

{Wl^a = + Pv9na + Pa9nv ■ (53) 

The only divergent terms here are C24, C35, and C36. 

We illustrate the Passarino-Veltman procedure by describing the evaluation of C\\ and C\ 2 in terms of Cq and the 
B functions. We start by multiplying Eq. (|52a|l by p± and p 2 : 

p\Cn +pi 2 C 12 = (q-pi)c = Ri , (54a) 
Pi2<7n +^2^12 = (q-p 2 )c = R 2 , (54b) 

where pij — pi -pj and ( )c is the integral operator on the RHS of Eq. (|5l"1) (so that for example Co = (l)c)- We write 
Eqs. (H2J as 

with the solution 



C12) \R 2 ) ' Vpi2 P2 
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We find R\ and R 2 by noting that 



where 



Then we see that 



Q-Pi = ~^{-{-q 2 + m 2 1 ) + {-{q+Pif + m 2 2 ) + h) , (57a) 
q-P2 = ^-{-{-{q + Pl) 2 + m 2 ) + {-{q + Pl +p 2 ) 2 + m 2 ) + f 2 ) , (57b) 



/1 = m 2 —m%+p\ , (58a) 
h = m 2 -m 2 + ( Pl +p 2 ) 2 -p 2 . (58b) 



Ri = (q -Pi)c = -^-{-B {m 2 ,m 3 ;p 2 ) + B (mi,TO 3 ;pi + p 2 ) 

+fiC Q (mi,m 2 ,m 3 ;p 1 ,p 2 )} , (59a) 

R2 = {q-Pi)c = -2-{--B (mi,m 3 ;pi +p 2 ) + B (mi,m 2 ;pi) 

+f2Co(mi,m 2 ,m 3 ;p 1 ,p 2 )} . (59b) 

Since the B functions are already known, only the scalar Co function remains to be computed. Similarly, the C 2 j 
functions can be evaluated in terms of the Cu's and -B's, etc. At each level in the ladder, only the scalar functions 
are new. 

The general three-point scalar integral is 

Co(m l7 m 2 ,m 3 ;pi,p 2 ) = J {dq)" [(~q 2 + m\)(-(q + pi) 2 + m 2 2 ) 

*(-(q+Pi+P2) 2 +^3)]^ 
dzdx^- , (60) 

where the limit n — » 4 has been taken since Co is ultraviolet finite, and (dq)" = (dq)'[ = d A q/(in 2 ). For A one finds 

A = (1 — z)m 2 + z(l — x)m 2 + zxm 2 — z(l — z)p\ — xz(l — z)2p\ 2 — xz(l — xz)p 2 . (61) 
The cases of interest here are 

C (0,m,m;puP2) = ^{M^^) 'M^)} ( 62 ) 

which holds when p\ = 0; 



C o (0,m,m;pi,p 2 ) = ^-{Li 2 (l - 2a) + 2C(2) - 2(arctan ^i— ^ ) 2 } 



(63) 



which holds when p\ — m 2 and (2pi + p 2 ) 2 = and where a = 2 + pi 2 /m 2 ; and 

CoK Wl ,p 2 ) = -±{L(k!i±££) - L(4)} (64) 
2pi 2 l V / V m / J 

where p| = and 

/" ln(l - z(l - z)s) / l~~s~\ 2 , , 

L(s) = y <fe v (1 z ) = -2(arctariW- J . (65) 

The dilogarithm function Li 2 (x) is discussed in detail by Lewin. |35| 
The four-point functions are defined by 



{D , Dp, D^D^vajD^ap} = J (dq)" {l,q t _ l ,q f _ t q lJ ,q tl q u q a ,q tl q u q a q !3 } 
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x [(-q 2 + mj)(-(q +pi) 2 + ml)(-{q+ Pl + p 2 ) 2 + m 2 ) 

x (~{q+Pi+p 2 +P3) 2 + ■ (66) 

We have dispensed with the regularization here since all of the D functions needed for our calculation are ultraviolet 
finite. While general expressions for Do exist, we need only a few special cases. In particular, we find that 

D (0,m,m,m;mN,-k 1 ,-k2)= . % = {Li 2 f-g±=, o) - Li 2 9) } , (67) 

where 

D = ^ , (68a) 
S12S23 

* ± = *(i±a/PP1 ■ (68b) 



S12S23 



tan( 




and Li2(r, 9) is the dilogarithm of complex argument. |35j | By some transformations among the momentum vectors, 
one can show that 

Z?o(0, m, to, to; mN — k±, — &2, — ^3) = Dq(0, to, to, to; mN, — £3, — hz) . (69) 

Finally, we also have 

D (m,ni,m,m;~k 1 ,~k 2 ,-k 3 ) = - _ 4 == (Li 2 (-7==, 0i) - Li 2 f — 1 

V S 12 S 12S23S23 L V V- ! 7 



%/Sl2Sl2S23S23 ^ "^VS7 / ~\y/Di 

+U 2 (^±=,9 3 ) -Li 2 (^,6 3 

\ V -L'3 ' ^ V i^3 

- Li2 (^'°) +Li2 (vft' )} ' (70) 

where 



D 


S12S23 

4S12S23 


; 


(71a) 


Di 


_ "1 2 S 23 
Sl 2 S 2 3 




(71b) 


D 3 


m 2 si 2 

Sl 2 S 2 3 




(71c) 


V± 


= i(i ±2. 


/Do") , 


(71d) 


0i 


= arctan * i 


I~S12 
S12 


(71e) 


0; 


— arctan * i 


/«23 
«23 


(71f) 



All of these integrals were done directly by way of Feynman parameters. 

The five-point functions are required for the ladder diagram. The five-point functions are very difficult to evaluate 
in general. We require only a special case, where mi = 0, m 2 = m% = m± = m$ = m, pi = mN, p 2 = —k\, p 3 = — k 2 , 
Pi = — k 3 . One feature of this special case is that there is a binding singularity: the scalar five-point function diverges, 
so we will have to base our implementation of the Passarino-Veltman formalism on the integral of the vector q^, 
which is finite, instead of on the divergent scalar integral. Also, we have not yet evaluated the three- and four-point 
functions with the necessary momenta. We give the three-, four-, and five-point functions with the special case mass 
and momenta values the names E, F, and G: 



(/)£ = ™ 2e j(€/[H 2 )(-(9+Pi) 2 +™ 2 )(-(9-Pi) 2 + ™ 2 )] 1 , (72a) 
</W) = ™ 2e f(dq):f[(-q 2 )(-(q+Pif+m 2 )(-(q- Pl f + m 2 ) 
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x (-(<? + Pi +p 2 ) 2 +m 2 )} 1 , (72b) 
</>G(p 2 , P3 ) = ™ 2 * j{d q )'U[{- q 2 ){-{ q +p 1 f + m 2 ){-{q-p 1 f + m 2 ) 

x(-(? + Pi +P2) 2 +m 2 )(-(g+pi +p 2 +p 3 ) 2 +m 2 )]~ 1 , (72c) 

where pi = p = mN. The first two of these are special cases of the three- and four-point functions. Because of the 
binding singularity, Eq — (1)e, Fq = (1)f, and Go = (1)g ah diverge. We start our analysis with the vector integrals 
Ef, = (q^E, = (q^p, and = (q^c- 

The three-point special case vector integral has the general form 

E» = (q») E =V»E 1 . (73) 

It is not hard to show (by use of symmetric integration) that Ei = 0, so that 

E ll = . (74) 

The four-point special case vector integral has the general form 

F„(p 2 ) = {q^F=Pi„F n (p 2 ) +p 2 „F 12 (p 2 ) . (75) 

The necessary vector integrals for p 2 — —k\ are 

Fn(-fci) = ^{h\ 2 {l-2x l )-2x l \n{2x 1 )- 26 2 -A^T 1 6 + 2C{2)} , (76a) 
F 12 (-h) = L Fll ^ kl ) + _L[2U 2 (l- 2x 1 )+((2)- 29 2 ] (76b) 



where 

When p 2 ~—k\ — k 2 one finds 
which implies that 



9 = arctan ^J— . (77) 

F^-h - k 2 ) = -F^-ks) , (78) 
Fni-h-h) = -F n (-k 3 ) + 2F 12 (-k 3 ) , 

F 12 (-k 2 -k 2 ) = F 12 (-k 3 ) . (79) 
The five-point special case vector integral has the general form 

G^-ku-kz) = (q„) G =P fl G 11 + k lf ,G 12 + h^G 13 . (80) 

The Gu functions are given by 

G11 (mi, a: 3 ) = -~^[lo{xi,x 3 ) + Ii(xi,x 3 )] - -^[I (x3,xi) +/i(x 3 ,xi)] , (81a) 
0X1 0X3 

Gi 2 (xi,x 3 ) = 1 _ [(1 - 2x 1 )I (x ll x 3 ) - Ji(xi,x 3 )] 
loxixi 

+ 16xxx 3 [ J o( a; 3^i) + h(x 3 , Xl )] , (81b) 
Gi 3 (xi,x 3 ) = 16 _ i 1 [lo{xi,x 3 ) + h(x 1 ,x 3 )] 

~ 2x 3)^o(^3^l) - M x 3,^l)] , (sic) 

/ (xi,x 3 ) = j _ [Li 2 (r+,fl)-Li 2 (r-,fl)] , (82a) 
VX1X1X3X3 

7i(xi,x 3 ) = rlnf— ) -^L= arctanfw — \ , (82b) 

(xi-x 3 ) Vx 3 / y/x 3 x 3 \y x 3 ) 



where 



(xi-x 3 ) \x 3 / ^/x 3 x 3 
with r± = ^fx~i ± \/ x\x 3 jx 3 and 9 — arctan \Jx\jx\ 
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:> 



\AAAA» 



\J »- 



\A/W\Jf- 







(a) 



(b) 



(c) 





(d) 



(e) 




(f) 



FIG. 2: Graphs contributing to the o-Ps decay amplitudes through order-a. They are the (a) self-energy, (b) outer vertex, (c) 
inner vertex, (d) double vertex, (e) ladder, and (f) annihilation contributions. The wave function factors are implicit in these 
graphs. 



VII. ANALYSIS OF THE ONE-LOOP DECAY DIAGRAMS 



The decay amplitudes can be written as 



Ai = A 



(o) 



A 



(i) 



A 



(2) 



(83) 



for i = l,2, 3, where the superscript indicates the power of a above that of the lowest order amplitudes = Af°. 

(Terms of order and higher also contain factors of ln(l/a).) The expressions for the squares \M\ lt \ 3t \ 3 - tm \ 2 contain 
parts of the form 



A*Aj 



[4 0) *A?> + A^*Af ] + A?>*Af> + [A? 1 * A 



i(i)* a(o) 



(i)* 4 (i) 



|(0)* A ( 2 ) 



Af^Af] 



(84) 



for various combinations of i and j. The A^*Aj terms give the lowest-order differential decay distribution. The 



, A^*A^ terms give the order-a correction, and the A) 1 '" AY' and AY 1 ' A, 

3 * J L J L J 

order-a 2 corrections. 

The graphs contributing to the order-a corrected decay amplitudes A^ in the renormalized Feynman gauge are 
shown in Fig. 2. The infrared divergence induced by mass-shell renormalization is regulated by use of a photon mass 
A. The self-energy (Fig. 2a) and vertex graphs (Figs. 2b, 2c) contain infrared divergences of the form lnAM^o- The 
ladder graph Fig. 2e requires special care in its evaluation since it contains an infrared binding singularity. This 
divergence can be identified and subtracted out, as discussed in detail in Ref. 0]. The result is that 



j 



|(0)* A {2) 



Ap^* Ay terms give the 



Mr = 



lnA-1 



Mil 



M. 



LO 



LS 



The subtracted ladder graph is 



M LS = 



with p = mN, 



and 



■ 4_2 

-ia m 



j (di)"[£ 2 (f - 2£p)(£ 2 + 2ip)Z{l)] 1 {{tv{l) - tr(0)) - f|j(^) - Z(0))} 

•S3 



tr(£) = -tr[^(j(£-p) + m)je^(£-p + k 3 ) + m)'ye* 2 (j{e + p~k 1 )+m) 
xjel(-f(£ + p)+ m) 7M ( 7 7V + l) 7 e( 7 iV - 1)] , 

Z{£) = ((£-p + k 3 ) 2 - m 2 )((£ + p - k^ 2 -m 2 ) . 



(85) 



(86) 



(87) 
(88) 
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The subtraction in Ea. l|86(l takes away the ^-independent part of ti(£)/Z(£), which would have had an infrared 
singularity. This binding singularity regulated by the photon mass, is displayed in Eq. JSSj. The contributions of the 
order-a decay graphs were evaluated one by one and summed. The 1/A binding singularity was removed according 
to the usual procedure of NRQED l36j . The In A terms cancel between the self-energy, vertex, and ladder graphs. 
The remaining expressions are a finite sums of rational functions of the Xi times logarithms, dilogarithms, and inverse 
tangent functions. 



VIII. RESULTS AND CONCLUSIONS 

We use our analytic results for the order-a decay amplitudes A± to calculate the order-a correction to the o-Ps 
37 decay rate and a part of the order-a 2 correction. The individual amplitudes are quite lengthy and will not be 
displayed. A simplified form for the complete order-a decay rate contribution is given in the Appendix. The result 
for the order-a decay rate is [32} 

Ti = -10.286606(10)-r LO . (89) 

This represents a 60-fold improvement in precision over the previous best result —10.2866(6) [3^ done using a higher 
dimensional integration. The two-dimensional integral for the part of the order-a 2 correction to the decay rate coming 
from the A^*Aj terms gives [39^ 

r 2 (square) = 28.860(2) (-) r LO • (90) 

The previous result for this contribution was 28.8(2). 

In this work we obtained analytic expressions for the o-Ps^ 37 decay amplitudes. We used these expressions to 
obtain precise results for the one-loop and "square" decay rate contributions, which were incorporated into the full 
calculation of two-loop corrections to the o-Ps— > 3j decay rate. [ToL 12^ We also give an explicit form for the one-loop 
decay distribution (see the Appendix) which can be used to obtain the one-loop energy spectrum in a convenient 
form. 
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APPENDIX: THE ONE-LOOP CORRECTION 



In the appendix we present the integral for the one-loop correction to the decay rate in compact form. From 
this integral the one-loop phase-space distribution and energy spectrum can be obtained. We note that for very 
soft photons additional effects must be taken into account in order to obtain accurate results for the phase-space 
distribution and energy spectrum. |4lJ, l42l l43l l44| 

The one-loop correction to the decay rate is 

ma 7 f 1 f 1 1 

ri = — „ / dx\ I dx2 \F(xx, x 3 ) + permutations) , (A.l) 

36?H J J 1 __ Xi X1X2X3 

where X1+X2+X3 = 2 and the "permutations" are the six permutations of the variables x%, X2, x 3 . The one-loop phase 
space distribution is just the integrand. The energy spectrum is found by integrating over X2 but not x\ — Ei/m. 
(The corresponding lowest-order expression is given in Eq. (1421) . The function F(x%,xs) is given by 



5 7 

F(xi,x 3 ) = go(xi,x 3 ) + } j gi(xi,x 3 )hj(xi) + / ] gi(xi,x 3 )hj(xi,x 3 ) 

z—1 i—6 



(A.2) 



The h functions are given by 

hi(xi) = ln(2xi) , 



h 2 (xi) = x j^-0i , 
V xi 

h3{Xl) = 2^{ C(2)_Li2(1_2xl) } 



h 5 (xi) 



1 P 



h 6 (x 1 ,x 3 ) = _ _ \lA 2 {rX,6i) -U 2 {r A ,6i)} 

VX1X1X3X3 I J 

h 7 (xi,x 3 ) = - = _ Ju 2 {r^,9 1 ) -Li 2 {r B ,e 1 ) 
^/xjxjxjxs I 

-±U 2 (r+,0) + hi 2 (r-,0)} , 



where Xi = 1 — Xi and 



9 1 = arctan(-\/ xi/xi 
0~i — arctan(-\/xi/xi 



V V hits/ 



V V ^1^3 / 

r|/v^r • 



.31 

.42 



The g functions are given in terms of x mn = x^x^ and x 2 = x\ + x 3 — 1 as 

^ 9x^(1- 2x0x3.3(1 - 2x 3 ) {-180 + 219^-4968^ 
+5292x 30 - 2 6 64x 40 + 504x 50 - 5848x n + 22639x 21 - 20280x 
+8405x 41 - 1240x 51 - 24x 61 - 175 51x 22 + 22 9 82x 32 - 5857x 4: 
+264x 52 + 48x 62 - 3776a; 33 - 878x 43 + 400x 53 + 536x 44 | , 

9i(xi,x 3 ) = -J7J ^— 7o7 — |2x 20 - 13x 30 + 35x 40 - 36x 50 + 8x 60 

xj(l - 2xi) 2 (xi - x 3 )x 3 I 

+4x 70 + 9a; 11 - 59x 21 + 149x 31 - 210x 41 + 162x 51 - 51x 61 + x n 
-4x 02 + 3x 12 + 55x 22 - 126x 32 + 104x 42 - 39x 52 + x 62 + 8x 03 
-26x 13 + 7x 23 + 22x 33 + 2x 43 - 2x 53 - 4x 04 + 14x 14 - 8x 24 - 8x 34 } 

92{xi,x 3 ) = —J (-48x 10 + 180x 20 - 2 76x 30 + 2 28x 40 - 108x 50 + 24x 60 

3xfxiX3 L 

+48x 01 - 48X 11 - 144x 31 + 244x 41 - 106x 51 + 2x 61 + 4x 71 - 96x 02 
+156x 12 - 108x 22 + 168x 32 - 132x 42 + 7x 52 + 6x 62 + 48x 03 - 60x 13 
-36x 23 + 42x 33 + 9x 43 - 6x 53 + 6x 34 - 4x 44 1 , 

93 ( Xl ,x 3 ) = 2/ 4 - |-2x 20 - 2x 40 - 4x 60 + 5X 11 - 6x 21 + 14x 31 - 4x 41 
xf(xi - x 3 )x 3 L 

+18x 51 - x 61 - 4x 02 - 2x 12 + 4x 22 - 2x 32 - 26x 42 - x 52 + 8x 03 
-7x 13 - 2x 23 + 12x 33 + 2x 43 - 4x 04 + 4x 14 } , 

gi {x u x 3 ) = A(_4 + 7x 10 - 7x 20 + 12x 30 - 10x 40 + 2x 50 + 8x 01 - lOx 11 
xf I 
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+3x 21 - 3x 31 + 2x 41 - 4a; 02 + 3x 12 + 2a; 22 + a; 32 } , (A.5e) 

g 6 (xi,x 3 ) = — \8-34x w + 29x 20 -4x 30 +6x 11 + 8x 02 -ix 12 ) , (A.5f) 
x\ V > 

96(xi,x 3 ) = — r 1 16 - 76a; 10 + 136a; 20 - 124a; 30 + 64a; 40 - 16a; 50 - 60a; 01 

X1X2X3 L 

+2723; 11 - 424x 21 + 294a; 31 - 104a; 41 + 22a; 51 + 92a; 02 - 392a; 12 

+484a; 22 - 187a; 32 + 13a; 42 + 2a; 52 - 76a; 03 + 294a; 13 - 259a; 23 + 30a; 33 

+3a; 43 + 36a; 04 - 120a; 14 + 61a; 24 + 3x 34 - 8x 05 + 22a; 15 + 2x 25 } , (A.5g) 

97(x 1 ,x 3 ) = ^-(l6-48a; 10 + 46a; 20 -12a; 30 - 2x 40 - 48a; 01 + 60a; 11 +9a; 21 
x 2 L 

-31a; 31 + 10a; 41 + 46a; 02 + 9a; 12 - 42a; 22 + 11a; 32 - 12a; 03 - 31a; 13 

+lla; 23 - 2x 04 + 10a; 14 } . (A.5h) 
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